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А. В. Красавин, А. Н. Чукарин, Б. М. Флек 
(Донской государственный технический университет) 


Несмотря на различия технологических процессов, реализуемых на фрезерных и сверлильных станках, пред- 
ставлен общий подход к построению моделей виброакустической динамики колебательной системы «шпин- 
дельный узел — режущий инструмент». Кроме этого на модельных фрезерных станках помимо традиционно- 
го фрезерования производится растачивание и сверление. Учитывая конструктивные особенности и геомет- 
рические параметры режущего инструмента, в качестве акустических моделей излучателя приняты линей- 
ный и точечный источники. Определение собственных частот и скоростей колебаний основано на использо- 
вании дифференциальных уравнений изгибных колебаний систем с распределёнными параметрами и функ- 
ЦИЙ А. Н. Крылова. Получены аналитические зависимости для скоростей колебаний, учитывающие парамет- 
ры технологического процесса, конструктивные особенности инструмента и шпиндельных узлов. 

Ключевые слова: виброакустическая динамика, шпиндельный узел, режущий инструмент, сверлильные и 
фрезерные деревообрабатывающие станки. 


Введение. Конструктивные особенности данного типа станков, заключающиеся в отсутствии 
зубчатых колёс в приводе, высокие частоты вращения шпинделей, многообразие типов фрез поз- 
воляют предположить, что в формировании спектров шума доминируют шпиндельная группа и 
режущий инструмент. Поэтому в качестве моделей акустических излучателей приняты: 

» линейный источник (при К > 344/1, ), уровни звукового давления которого на основании 


данных [1, 2] для рассматриваемых станков приведены к следующему виду: 
Е =20 щи, +2081, + 2019/ +109 А —- 2099г +124; 
» точечный источник (при А < 344/Ё, ), звуковое давление которого приведено к следую- 


щему виду 
Е =201щи, +2082, +2014! - 2014г +106, 
где и, — скорость колебаний источника; & — собственные частоты колебаний источника; А — ра- 
диус шпинделя или фрезы; / — длина источника; г — расстояние от источника до расчётной 
точки. 
Результаты исследований. Как видно из полученных зависимостей для расчёта уровней звуко- 
вого давления следует определить собственные частоты колебаний и виброскорости на этих ча- 
стотах. Для этого в работе, в отличие от существующих, рассмотрены три варианта, приведённые 


на рис. 1. 
| | | | 
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Рис. 1. Расчётная схема определения виброскоростей системы «шпиндельный узел — инструмент» деревообрабатываю- 
щих станков: а — отрезание или обработка торца изделия; 6 — фрезерование паза; в — схема, учитывающая весь шпин- 
дельный узел; г — упрощённая схема фрез малого диаметра 


1 Работа выполнена в рамках инициативной НИР. 
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Определение собственных частот колебаний и скоростей колебаний основано на исполь- 
зовании дифференцированных уравнений изгибных колебаний систем с распределёнными пара- 
метрами и функций А. Н. Крылова. 

Общее дифференциальное уравнение изгибных колебаний имеет вид: 

д д22 078 
В, т 
Хх 
(1) 
д 5 д27 
ОВС = 
ох ох ОЕ 
где Е — модуль упругости, Па; 7 — момент инерции в направлении соответствующей оси, м“; 
т — распределённая масса, кг/м. 

Решение уравнения, соответствующее характеру возбуждающей силы, которая фактиче- 

ски является функцией периодического характера, определяется из уравнения: 


92(х) 4 
——^-а 2(х}=0 
дх4 (к) 
д%у(х) 4 
———^-а ух) =0 
дх4 (к) 
тождественного выражению свободных колебаний при значении 


2 РИ п 
а = — ти 
Г 1, 


где Г. — собственные частоты колебаний, Гц; № — частота силового воздействия, Гц. 
Первая схема. Решения уравнения (1) на участках фрезы, свободных от нагрузки, представим в 
виде [3]: 
2(х) = АК, (ах) + ВК. (ах) + СК-(ах) + Ок. (ах), (2) 

1 . 1 . 

где А, (ах) = 1 (спах + созах}; К› (ах) = 1 (пах + пах); к, (ах) = 1 (спах + совах ); 
Вы 
Ка (ах) = (пах + пах} 

Постоянные интегрирования определим из граничных условий. В точке приложения 
нагрузки должны выполняться условия сопряжения, требующие равенства прогибов, углов пово- 
рота и изгибающих моментов для обоих участков фрезы, а также скачка поперечной силы, равно- 
го по величине возмущающей силе. Поэтому 

Р 
2(х) =, =2 Иа + Е: К. [а (х- /)]. 
Располагая начало координат на правой опоре, получим 
2 (Х) ов = ВК, (ах) + ОК, (ах), 


так как постоянные Аи С равны нулю в соответствии с граничными условиями на правой опоре. 
Соответственно на левом участке амплитудные значения смещений выражаются формулой 





2(х)лев = ВК. (ах) + ОК. (ах) + в ках -Ь)]=0. 





а3Е] 
Постоянные Ви ДР определим из граничных условий на левой опоре 
а2(х) _ р, 
о ВК(а/)+ ОК. (а!) + Е ках -Ь)= 0, (3) 
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. а ви р (©/) + 55 Кох > О 





4?2(х) 
ах? 
Из системы (3) определяем 


в-_Р |“ = 5) (о!) ® [о - 5 )№(@/) Ка (@/) _ а [6 и 








аЗЕ7 К2(а/)- К2(@/) `ю(@/) — Ю(@)) 
р___Р Ка - 5) (©!) - [а — 5 (а!) 
аЗЕ) кз (а/)- К? (@/) | 


Тогда величина прогиба по оси О2 равна 
2(х,Ё) = [ВК. (ах) + ОК. (ах (2). 
Аналогичным образом находится прогиб и по оси ОУ. 
02 ду 


Радиальная скорость колебаний определяется формулой —. и т 


Выражения для у(х,Ё) и бу (к,Ё) имеют аналогичный вид с учётом замены Р, и Р, 


ОЕ 


на в и Ра 


Собственные частоты колебаний заготовки определяются по известной формуле 


22 0,5 
ды "|=. ("=123,..) 
2/^ То 
Полученные выше зависимости справедливы для абсолютно жёстких опор. Фактически 
передний подшипник шпинделя с патроном и обрабатываемая заготовка являются упруго- 
диссипативными опорами, податливость которых существенно влияет на динамические парамет- 
ры системы «шпиндель — заготовка». 
Вторая схема. Рассмотрим теперь фрезу, опёртую по концам на две упругие опоры с коэффици- 
ентами податливости с, (заднего центра) и © (переднего шпиндельного подшипника и зажимного 
устройства). Начало координат расположим на заднем центре. Граничные условия на концах за- 
готовки запишем в виде: 
х=0 у= Ос; а?у/ах = 0 7 
х=!/ у= 06; а?у/ах = 0, (9) 
где О — амплитудное значение поперечной силы в соответствующем сечении, определяемой по 
формуле 
4°у 
ахз 
С учётом последнего выражения получим граничные условия в следующем виде 


азу 
| + СЕ] а 


9=Е 


=0 





х=0 
у 


ах? =. 


х=0 


аЗу 
| с›ЕЛ 3 








=0 





Ж= 
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Подставляя в эти уравнения общее выражение для прогибов (2) и учитывая соотноше- 
ния функций А. Н. Крылова, получим следующую систему уравнений относительно постоянных 


АВ, СР 
А+сЕХЗЬ =0, 
К, (^) - с»ЕЗазк, (^) |+Р [К (^)- с.Елазк, (Л) | =0, 


А[ К, (^) с,ЕЗабК, (Л) | + В[ К, 
АКз(^) + ВК. (Л) ОЖ (Л) = 0, 


где Л =а.. 
Решая данную систему, получим выражение для определения собственных частот 





колебаний 
($НЛ — Л). — с,Е?а? (СНА + созл) ($НЛ — Л) — СИЛ — с0$А — с,Е?а? (эНЛ + Л) с,ЕЛа? ($ВЛ — пл) 
— [$АА + $1 — с,ЕЛа? (сВЛ + созА) || $ВЛ + эл — с,Е?а? (СПА — соз\) | = 0. 
Представляя нагрузку с помощью дельта-функции и учитывая краевые условия (4), полу- 
чим дифференциальное уравнение поперечных колебаний заготовки на податливых опорах 
Зплх ] 
=" — 











о ^ 1 пих 
+а = Р.УЕ Ах +В +—| 35 
0? дхА ти у РЗ 4 / 
С.В, - С.К С.К 
где А=-? г 11, = Г т, К, и Е, — реакции в опорах заднего центра и в переднем шпин- 
2 
ё Р, и Р, на 


дельном подшипнике. 
Аналогичным уравнением описываются колебания по оси у с учетом замены 


Рир. 
Эти уравнения в конечном виде не интегрируются, поэтому для их решения использованы 


численные методы. 
В качестве начальных условий принимаем (х =/ и Ё =0) 


у=с,Р, и ау [а =а(с В, /а&, 


где К› — реакция в переднем шпиндельном подшипнике 
Решение этой задачи с использованием функций А. Н. Крылова имеет вид 


Р 
Улев = -Е С.К, (ах ) + ВК, (ах) + ОК. (ах) + ЕТ [а(/- 5.) |. 


Постоянные Ви р определим, исходя из краевых условий на левом конце заготовки, получаем 


К. (а/), 


систему уравнений: 
Р 
-=Ю с. =-ЮСК (а/) + ВК, (а/) + Ок, (а/) + РЕ 




















Я 
ЧУ -0- Ас, (9/)+ 86, (07) + 0%, (0) + 5, (61), 
из которой вычисляем 
Таотиа „5 о И ‚ (а/) Р К а(!-51) К, (а/) + К, [а(/- 54) К, (а!) 
О= | ‚ 
К, (а!) К: (а/) т м ЕЛ К; (а/) + К (а!) 
р, К а( __ к, (а/) + к, [а(/ - 58) | К, (а/) К, (а!) Р,к, [а(1- 58) | 
8-Е. Ко)? (а/) 1. (67) — <ЕЖ, (67) 
. ) + К,с, ‚ети К, (а/) вс (4!) 
Ка (а/) + К? (а!) "ТК, (а!) 





№ (1) (а!) 
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Вибросмещения по оси ув этом случае определяются выражением 
ух, Е) = вс. (ах)+ ВК (ах)+ Ока (ах) + С°япоЕ) 
Теперь определим радиальную скорость колебаний как >. Аналогичными выражениями описы- 


92(х,Ё) т Е 


ваются Хх, В и с учётом замены Р, и Рна Р, ир. 


В окончательном виде прогибы режущего инструмента по оси г определяются выражением 


2 (х) = [А (спах + созах ) + В (пах + пах ) +20 (пах _ пах ) | ПЕ. 


Прогибы режущего инструмента по оси ОУ находятся по аналогичному выражению, в ко- 
тором сила Р‚ заменяется на Р.. 
Третья схема. Для существенного уточнения колебательной модели следует рассмотреть систе- 
му шпиндель — заготовка, которая представляет собой балку на двух опорах с консольной ча- 
стью. В данной расчётной схеме приняты следующие обозначения: — средневзвешенный мо- 
мент инерции межопорной части шпинделя, а № — средневзвешенный момент инерции консоль- 
ной части. 
Для каждого из участков запишем выражения прогибов, используя функции А. Н. Крылова: 
у! = «Аа (Ах) с>К5(Ах)+ сзКз(Ах)+ саКа(Ах) при 0 <х <, 
уз = с! (Ах) + ск (Ах) сзКЗ(Ах) + с.Ка (Ах) при < х</, 
где /— длина шпинделя и заготовки, м; Ь — длина межопорной части шпинделя, м. 
На первом участке у, (6) = у. (0)=0 
Поэтому с! =с; =0 
В точке сопряжения двух участков 
уз (6) = у), т 9 | 
х=В х=В 
Поперечные силы в конце первого и начале второго участка отличаются на величину опорной 
реакции. 
Поэтому 








93 2 
Е Е? К. (5) 
ох 9х? | _ 
х=Ь х=Ь 
Условия сопряжения участков можно выполнить, если представить прогибы на втором 
участке в виде 














В 
72у> = ЛУ; + ЗЕ Ка [Ах ь) (6) 
Функция Ал в выражении (6) тождественно равна 0 при х<Б и не равна 0 при х>Б. Послед- 
нее выражение справедливо для всей системы «шпиндель — заготовка». Подставив в (5) выра- 
жение для у получим (с = < = 0) 


У К 
= 5. Ех) ска (Ах +3 ах -6)] 


0,5 
В 
пе ^=[*% | р 


Полученное выражение для у содержит постоянные ©, с! и К, которые можно определить 
из условий 
Х=Ь у =0 
х=/ уУ=Уу=0. 
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В итоге получаем систему из трёх уравнений 
ско (АЬ) + са (АЬ) — 0, 
Г 


Е № [Аи -Ь)]=0, (7) 





=. [Ка (^№)+ са (+ 
2 


ом) см БМ 6] =0. 


Приравнивая нулю определитель системы (7), получим уравнение для расчёта собствен- 
ных частот колебаний: 











ИЛЬ + УпЛЬ ИЛЬ — эпЛЬ 0 
21 (1 — яп) 2 (вм ам) АР) ее 
Л, 7, А Е, 
зам азам: РОВ 
| 7, А Е, | 





Четвёртая схема. Координата приложения нагрузки на фрезе не меняется относительно точки 
закрепления. Расчётной моделью является консольно-защемлённая балка. Прогибы фрезы по 
осям О2и ОУопределяются из уравнения (1), в котором 


2 2 
2 @П [р 
а = —— . 
р 215 





Здесь #„, — собственные частоты колебаний фрезы, Гц; /, — длина консольной части, м. 
Для функции 2(х) = АК, (ах) + ВК. (ах )+ СК (ах + Ок. (ах) имеем граничные условия 
Р.[3 
РНЕ: ‚ 2 
ЗЕЛ 





При 2(х =0)=0 получаем, что С= 0, а 


2(х =0)=А= 





Из остальных граничных условий находим: 
— АК, (а/ + ВК>(а/›)+ Ок. (@/,)= 
— АКз(а/›)+ ВК. (а/,)+ ОЖ (а/,}= 
После преобразований получаем 
_ РЫр спа/ьпа/» — зва/,соза/, 
6Е7 р эпа/,5та/, 


0 
0 








Й 


_ РУз ста/ па!» + эва/соза/, 

6Е7„ эпа/,5та/, 
Выводы. Приведённые результаты исследований позволяют существенно уточнить модели 
виброакустической динамики основной колебательной системы «шпиндельный узел — режущий 
инструмент» модельных и сверлильных деревообрабатывающих станков. Обоснован общий под- 
ход к расчёту процесса шумообразования этих станков. Различия в самом расчёте заключаются 
только в задании силового воздействия от процесса резания. Полученные зависимости учитывают 


В 
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все основные конструктивные параметры как режущего инструмента, так и шпинделя, а также 
параметры технологического процесса. 
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;триеа зуяетс апа А. М. Ктуоу$ ипсВоп5. ТВе апа&с аерепаепсез Юг те ифгайоп уеосШез сопеттд, рго- 
се55 уапаЬ/Ез, дез!ап Геайигез оЁ {те юо/ апа рт ипй5 аге оматед. 

Кеуигога&: ибгоасоиз Вс аупати/с$, рта ит тас@тта юо, тИтд апа аита иоодИ/огКЕГе. 


1 Тпе гезеагсИ {5 допе мййт {пе Я’ате о! {Не паерепдеге В&0. 
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